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Abstract. The protostellar outflow HH211 is of considerable interest since it is extremely young and highly collimated. Here, 
we explore the outflow through imaging and spectroscopy in the near-infrared to determine if there are further diagnostic 
signatures of youth. We confirm the detection of a near-infrared continuum of unknown origin. We propose that it is emitted 
by the driving millimeter source, escapes the core through tunnels, and illuminates features aligning the outflow. Narrow-band 
flux measurements of these features contain an unusually large amount of continuum emission. [Fell] emission at 1.644 /an 
has been detected and is restricted to isolated condensations. However, the characteristics of vibrational excitation of molecular 
hydrogen resemble those of older molecular outflows. We attempt to model the ordered structure of the western outflow as a 
series of shocks, finding that bow shocks with J-type dissociative apices and C-type flanks are consistent. Moreover, essentially 
the same conditions are predicted for all three bows except for a systematic reduction in speed and density with distance from 
the driving source. We find increased K-band extinctions in the bright regions as high as 2.9 magnitudes and suggest that the 
bow shocks become visible where the outflow impacts on dense clumps of cloud material. We propose that the outflow is carved 
out by episodes of pulsating jets. The jets, driven by central explosive events, are responsible for excavating a central tunnel 
through which radiation temporarily penetrates. 

Key words. ISM: jets and outflows, stars: circumstellar matter, infrared: ISM, ISM: Herbig-Haro objects 



1. Introduction 

The processes controlling the birth of individual stars evade di- 
rect analysis in the optical because they operate within highly 
obscuring cores within dense molecular clouds. Nevertheless, 
the end of the period of gestation is heralded by often spectac- 
ular ejections of material JBachillerJ fl996). We may thus de- 
duce information concerning the accretion process by indirect 
means. Bipolar jets from protostars appear to be ubiquitous and 
their characteristics suggest a close rel ationship to the accretio n 
process from a disc into a forming star jFroebrich et alJ20*0 3b). 
In addition, they have drastic dynamical and chemical effects 
on their environment. They excavate cavities by sweeping up 
the molecular material which then forms an outflow visible in 
emission from CO molecules. As well as a leading bow wave 
which ploughs through the ambient medium, irregularities in 
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(CSIC). Observations were also obtained with The United Kingdom 
Infrared Telescope, which is operated by the Joint Astronomy Centre 
on behalf of the U.K. Particle Physics and Astronomy Research 
Council. 



the outflow give rise to shocks which propag ate away f rom 
the pro tostar into the ambient molecular gas dFroebrich et alJ 
2003a). These bow shocks excite atoms and molecules, in- 
ducing line emission which is often visible at optical wave- 
lengths as Herbig-Haro (HH) objects jReipurth & Balrvll200lt 
Bally & Reipurth However, the near-infrared emission 

is particularly important for the youngest protostars since even 
the protrud ing outflows may s till be obscured in the opti- 
cal regime (Eisloffel et al. 200$. In fact, near-infrared obser- 
vations reveal remarkable molecular hydrogen flows which 
hold valuable information concerning the environment through 
which they prop agate as well as the mechanism which launches 
the driving jets (Davis 2002). 

HH 211 is a bipolar molecula r outflow which was discov- 
ered bv lMcCaughrean et alJ (fl994). It lies near the young stel- 
lar cluster IC 348 IR in the Perseus dark c loud complex at an 
estimated distance of 315 pc (Herbig 1998). The outflow is bi- 
laterally symmetric and highly collimated with an aspect ra- 
tio of ~ 15:1. The total extent of the outflow is 106" which 
is 0.16pc at the adopted distance. A H2 (1,0) S(l) wide-f ield 
survey of the IC 348 cluster was earned out by Eisloffel et al. 
(2003) covering a 6'. 8 x 6'. 8 region; no HH 211 outflow rem- 
nants were detected beyond the outer knots. This makes it one 



2 



O'Connell et al.: The HH 21 1 protostellar outflow 



of the smaller outflows, which suggests that it may also be one 
of the youngest since the average length within an unbiased 
samp le of Class O/I jets was found to be 0.6-0.8 pc JStankd 
2003). 



The conclusion that HH 211 is a jet driven outflow with a 
timescale of order 1000 years was deriv ed from interferometric 
CO observations dGueth & Guilloteaull 19991) . which confirms 
it as one of the youngest infrared outflows to be discovered. 
Although the outflow lobes are visible in both blue-shifted and 
red-shifted CO emission, a small inclination angle to the plane 
of the sky is suggested by (i) the lack of strong differential ex- 
tinction in the H2 brightness distribution, (ii) the high degree 
of separation of the blue and red CO lobes and (iii) the rela- 
tively small radial components of the SiO a nd H2 flow speeds 
dChandler & Richerl200ltlSalas et all2003l 

The central engine driving the outflow is HH 211-mm, a 
low-mass protostar with a bolometric luminosity of 3.6 Lq and 
bolometric temper ature of 33 K. It is surrounded by a ~ 0.8 Mq 
dust condensation (Froeb richl20"05l) 1 . Since 4.6 per cent of the 
bolometric luminosity is attributed to the submillimeter lumi- 
nosity, L smm , HH 21 1-mm is classified as a Class O type proto- 
star. 

A compact and collimated SiO jet extends in both flow 
directions o ut to a projected distance of 20 " from the cen- 
tral source dChandler & RicheJll997l l200lb . This is also an- 
other indication of the Class O nature of the source JGibb et alJ 
120041) . The clumpy nature of the observed SiO emission sug- 
gests a shock origin re sulting from a time dependent jet veloc- 
itv. lNisini et alJ (120021) observed HH 21 1 in SiO lines originat- 
ing from high rotational energy levels and deduced a high jet 
density of n#2 ~ 2-5 x 10 6 crrT 3 and gas temperature > 250 K. 
However, no SiO emission is detected beyond 20"where the 
H2 (1,0) S(l) emission is found suggesting that the conditions, 
such as shock velocity and pre-shock density, vary consider- 
ably along the flow direction. 

Imaging of [Fell] emission in the H-band at 1 .644yum has 
recently playe d a major role i n our understanding of shocked 
outflows (see iReipurfh et ail feOOOh for a summary). Where 
[Fell] emission is observed it traces the fast (> 50 km s _1 ) 
and dissociative shocks. Combined with K-band imaging of 
H2 rotational-vibrational lines, which trace less extreme shock 
conditions, import ant information about outflow excitation 
may be gathered ( Khanzadv an et alJ2004h . 

Previous studies employing bow-shock models have fo- 
cused on larger, more evolved sy stems (Eisloffel et al. 2000; 
ISmith et all2003tlO'Connell et alJ2004 . HH 21 1 provides us 
with the unique opportunity of studying what has been deemed 
as an exceptionally young outflow covering a small spacial ex- 
tent. We present here new high resolution images of H2 and 
[Fell] lines and K-band spectroscopy of the outflow (§ 3). We 
analyse the impact regions visible in the near-infrared. We then 
present in § 4 C-type bow-shock models in order to interpret the 
remarkable set of bows propagating through a changing envi- 
ronment along the western outflow. In § 5 we discuss the issues 
which have arisen from this set of data before pooling our find- 



ings together in conjunction with previous studies to suggest a 
global outflow mechanism for HH 21 1 in § 6. 



2. Observations and data reduction 

2.1. KSPEC observations 

Our near-infrared (NIR) spectra cover the 1-2.5/^m region in 
medium resolution. They were obtained in the period 26-29 
August 1996 with the KSPEC spectrograph on the University 
of Hawaii 2.2 meter telescope. This cross dispersed Echelle 
spectrograph is equipped with a HAWAII 1024 x 1024 detector 
array and optimised for 2.2yum. Observations were performed 
at two bright H2 emission locations. The 0."96 width slits ran 
in an east- west direction passing through knots / and d in the 
west and through knots i and j in the east (positions are in- 
dicated in Fig.Q. Data reduction (flatfielding, sky-subtraction, 
extraction of the spectra) was performed using our own MIDAS 
routines. An absolute calibration of the fluxes was not possible 
due to non-photometric weather conditions. Wavelength cali- 
bration w as performed us i ng OH -night-sky emission lines and 
tables of iRousselot et alJ J2OOO). The results are presented in 
Tabled Note that fainter emission lines are detected in knots i 
and j due to the stronger emission here compared to knot d (see 
Table|2l) 

Table 1. KSPEC relative fluxes. For both slit positions, the 
measured fluxes lie above the continuum and are presented rel- 
ative to the H2 (1,0) S(l) line flux. Three observations were 
carried out at each slit location. The relative fluxes have been 
averaged and the spread in values is quoted in brackets as an 
error estimate. 



Line 




east 


west 


(1,0) S(9) 


1.687 


0.02 (0.01) 




(1,0) S(7) 


1.748 


0.12(0.01) 




(1,0) S(6) 


1.788 


0.07 (0.01) 




(1,0) S(5) 


1.835 


0.62 (0.04) 




(1,0) S(4) 


1.891 


0.20 (0.01) 




(1,0) S(2) 


2.033 


0.36 (0.01) 


0.30 (0.05) 


(3,2) S(5) 


2.065 


0.03 (0.01) 




(2,1) S(3) 


2.073 


0.10(0.01) 


0.18(0.07) 


(1,0) S(l) 


2.121 


1.00 (0.05) 


1.00 (0.05) 


(3,2) S(4) 


2.127 


0.02 (0.01) 




(2,1) S(2) 


2.154 


0.05 (0.01) 




(3,2) S(3) 


2.201 


0.03 (0.01) 




(1,0) S(0) 


2.223 


0.25 (0.01) 


0.31 (0.02) 


(2,1) S(l) 


2.247 


0.13(0.01) 


0.17(0.06) 


(2,1) S(0) 


2.355 


0.03 (0.01) 




(3,2) S(l) 


2.386 


0.03 (0.01) 




(1,0) Q(l) 


2.406 


0.96 (0.05) 


1.25 (0.17) 


(1,0) Q(2) 


2.413 


0.38 (0.05) 


0.54 (0.11) 


(1,0) Q(3) 


2.423 


0.95 (0.05) 


1.13 (0.05) 


(1,0) Q(4) 


2.437 


0.33 (0.04) 


0.41 (0.15) 


(1,0) Q(6) 


2.475 


0.19(0.06) 




(1,0) Q(7) 


2.499 


0.43 (0.03) 





1 http://www.dias.ie/protostars/ 
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Fig. 1. Broad-band K image of H H 21 1 which cove rs wavelengths between 2.03//m and 2.37/mi. Features are labeled according 
to the nomenclature of iMcCauehrean et al] Jl994 . The exciting source HH211-mm at position R.A.(2000) = 03 h 43 m 56.7 s , 
Dec(2000) = +32° 00' 50.3" ( Avila et al. 2001 ) is indicated by the cross. KSPEC spectroscopic slit positions are indicated by the 
dotted lines. Note that the star in the south-west of this image (and following images) has been masked for display purposes. 



2.2. MAGIC observations 



The NIR images were taken in November 1995 at the 3.5 me- 
ter t elescope on Calar Alto using the MAGIC infrared cam- 
era llHerbst et alJI 19931) in its high resolution mode (0.32" per 
pixel). Images were obtained using narrow band filters cen- 
tered on the H2 (1,0) S(l) emission line at 2. 122 /mi, the 
(2,1) S(l) line at 2.248 //m, the (3,2) S(3) line at 2.201 //m 
and on the nearby continuum at 2. 140 /mi. The per pixel in- 
tegration time was 1740s. Seeing throughout the observations 
was ~0.9"except for the 2.14 /mi image where it is ~1".8. The 
narrow -band image containing the ( 1,0) S (l) line has been pre- 
viously published in Eisl offelet all J2003I) . 

The data could not be accurately flux calibrated due to 
non-photometric conditions. The total H2 (1,0) S(l) flux from 
the entire outflow was previou s ly me asured to be 1.0 x 10~ 15 
W m -2 by iMcCaughrean et al] i 19941) which agrees with our 
flux calibration. We have deduced a calibration for all our nar- 
row band images according to this measurement assuming that 
the average integrated flux from several bright unsaturated field 
of view stars should be similar for each filter. This should be 
the case because each filter FWHM is equal to 0.02 /mi and the 
spectral energy distribution (SED) is likely to be relatively flat 
on average between 2. 122 /mi and 2.248 /mi. We assume the 
accuracy of this method to be ~ 15%. 



2.3. UFTI observations 

Further NIR observations of HH211 were carried out on 
December 12, 2000 (UT) at the U.K. Infrared Telescope 
(UKIRT) using th e near-infrared Fast Track Imager UFTI 
(Roche et alj l2003). The camera is equipped with a Rockwell 
Hawaii 1024 x 1024 HgCdTe array which has a plate scale of 
0".091 per pixel and provides a total field of view of 92". 9 x 
92".9. 

Images in the [Fell] 4 D7/2 - 4 Fi)/2 transition were obtained 
using a narrow-band filter centered on A = 1.644 /mi with 
A/l(FWHM) = 0.016/mi). The outflow was also imaged us- 
ing the broad-band K[98] filter centered on A — 2.20/mi with 
A/l(FWHM) = 0.34/mi. Seeing throughout the observations 
was ~0".8. Nine-point 'jittered' mosaics were obtained in each 
filter. 

Standard reduction techniques were employed (using the 
Starlink packages CCDPACK and KAPPA) including bad- 
pixel masking, sky subtraction and flat-field creation (from the 
jittered source frames themselves). The images were registered 
using common stars in overlapping regions and mosaicked. The 
observations were conducted under photometric conditions, so 
the faint standard FS 1 1 (spectral type A3; H = 11.267 mag) 
(Hawarden et al. 2001) was also observed and used to flux cal- 
ibrate the [Fell] image. The final images were smoothed using a 
circular Gaussian filter of FWHM = 0".25. in order to increase 
the signal to noise without compromising the resolution. 
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Fig. 2. HH21 1 at 2.248 fim. The grey-scale is the continuum subtracted image showing only the (2,1) S(l) line emission. Overlaid 
are contours representing the non continuum subtracted image. Thus the contours nearer the source trace predominantly scat- 
tered light while the contours further out mostly trace compact, line-emission features. The contour levels, which are scaled 
logarithmically, are at 0.22, 0.44, 0.87, 1.74, 3.47 x 10~ 18 W nr 2 arcsec" 2 . 
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Fig. 3. HH211 at 1.644 /im. The greyscale and contours both trace the [Fell] line plus continuum emission. Consequently, 
most of the distributed flux can be attributed to continuum emission, although some concentrated condensations of [Fell] line 
emission are observed. These are labelled 1 -3. The logarithmic contours measure the flux at 0.63, 1.0, 1.58, 2.51, 3.98, 6.31 x 
10" 19 Wnr 2 arcsec" 2 . 
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Table 2. Photometric results for HH 21 1. The flux in units of 10~ 18 WirT 2 is measured over the indicated circular apertures. Note 
that the measurements were made from images which are not continuum subtracted because of the low S/N of the continuum 
image. To infer the line emission fluxes the continuum values* need to be subtracted. The (2,1) S(l)/(1,0) S(l) (labeled 2/1) ratios 
have been determined after subtracting the 2.14yum continuum emission. In the case of knot g the (2,1) S(l) flux is equal to the 
continuum flux so no ratio was derived. 



Knot 


aperture 


(1,0) S(l) 


/"> 1 \ C / 1 \ t 

(2,1) S(l)' 


(3,2) S(3) 


O 1 A . f .. , t* 

2.14yum Cont. 


Broad-band KJ 


rT?~TT1 1 A A 

[fell J 1.644yum' 


2/1 ratio 


HH211-west 


57" 


413 


61 


34 


42 


1285 


28 


0.05 (0.04) 


HH211-east 


42" 


575 


98 


63 


55 


2197 


41 


0.08 (0.03) 


bow-a 


5.3" 


7.2 


0.6 


no det. 


no det. 


9 


no det. 


0.08 (0.04) 


bow-bc 


7.7" 


32.0 


2.3 


no det. 


no det. 


65 


no det. 


0.07 (0.02) 


bow-de 


10.1" 


155.5 


18.1 


7.3 


3.1 


335 


11.8 


0.09 (0.02) 


b 


4.9" 


19.3 


1.5 


no det. 


no det. 


41 


no det. 


0.08 (0.02) 


c 


4.9" 


11.2 


0.9 


no det. 


no det. 


21 


no det. 


0.08 (0.03) 


d 


6.0" 


89.6 


10.2 


4.0 


1.9 


192 


6.1 


0.09 (0.03) 


e 


6.0" 


57.5 


7.9 


3.3 


1.3 


131 


4.8 


0.12(0.03) 


f 


10.1" 


157.2 


19.3 


9.1 


8.9 


395 


9.0 


0.07 (0.04) 


g 


10.1" 


41.9 


20.6 


20.6 


19.6 


380 


14.3 




h 


10.1" 


60.3 


9.8 


5.1 


4.7 


190 


4.6 


0.09 (0.08) 


ij 


17.9" 


458.3 


54.3 


27.9 


24.6 


1364 


25.0 


0.07 (0.03) 



">' The background lcr noise estimates over an 8"aperture are: 2.8 for H 2 (1,0) S(l); 0.2 for H 2 (2,1) S(l); 0.7 for H 2 (3,2) S(3); 1.1 for the 
2.14/jm continuum; and 3.4 for the [Fell] at 1.644/jm, also in units of 10~ 18 WnT 2 . 

* As a standard star was not observed in the K[98] filter the broad-band fluxes are rough estimates. The flux calibration factor which we used 
was derived from the (1,0) S(l) calibration and a comparison of the broad-band and narrow-band filters. 



3. Results 

Fig- Q] displays the HH 211 outflow in the K-band between 
2.03/mi and 2.37 yum which contains all the K-band line emis- 
sion as well as a large proportion of continuum emission. The 



rincip le knots have been labeled as in McCaughrean et all 



1199 

The H2 (2,1) S(l) image is displayed in Fig. [2 The contin- 
uum at 2.14/vm has been subtracted in order to indicate loca- 
tions of pure (2,1) S(l) emission (greyscale). Contours of the 
non continuum subtracted image are also displayed in order 
to indicate the extent of the continuum emission at 2.248 /im. 
The line emission is produced from an excitation level which is 
12,553 K above the ground state whereas the (1,0) S(l) arises 
from 6,953 K. Therefore, we expect it to highlight the hotter 
parts of molecular shocks. 

The western outflow shows particularly interesting struc- 
tures which can be described as a series of bow shocks propa- 
gating along the outflow away from the source. Bows de and be 
display a common asymmetric structure: the lower bow wing 
is approximately 1.5 times brighter than the upper wing. We 
will interpret this asymmetry in §0]as due to a misalignment of 
the magnetic field with the flow through which the bow shock 
configurations with C-type flanks are propagating. 

Emission detected at 1.644yum is displayed in Fig. [3] Most 
of the emission detected here is attributed to the high level of 
continuum flux, as is seen in the K-band (Fig. [Q. Steeply ris- 
ing above this continuum level are several concentrated [Fell] 
emission condensations, labeled 1-3, one of which forms part 
of a well defined bow-shock, bow-t/e. 

We confirm the existence of a band of continuum emission 
which extends along the western outflow. It becomes visible 



5" from the driving source and ma intains a relatively con stant 
flux out to 17" from the source. Eisloffe Tet alJ d2003l) sug- 
gest that this continuum is scattered radiation from HH211- 
mm which opens the possibility of indirectly obtaining a spec- 
trum of the outflow source. The band of continuum is promi- 
nent in the K-band as well as at 1 .644/im. 

The photometric results for knots a - j are listed in Table|2] 
Knots d,e and b,c are interpreted as bow-shock components and 
are labeled as bow-ate and bow-bc. We will briefly discuss the 
implications of these fluxes. The eastern outflow is ~1.5 times 
brighter than the western outflow. However, the cause of this 
di fference is not necessa rily due to an unequal jet power output 
as lGiannini et alJ J200ll) detected similar levels of Ol 63 jim in 
both lobes (1.02 Lq in east; 0.94 Lq in west). This line is rel- 
atively unaffected by extinction and represents the main cool- 
ing channel in the post-shocked gas. The non photometric con- 
ditions during the observations are reflected in the relatively 
large lcr noise levels which have been estimated from the back- 
ground flux variation using an 8"aperture. Most of the emis- 
sion in the (3,2) S(3) image is actually continuum emission ex- 
cept possibly for bow-c/e. Interestingly, the fluxes at 1.644 yum 
are also comparable to the continuum fluxes. This implies that 
the dereddened continuum fluxes are 2-5 times brighter at 
1.644pm than at 2.14/im taking into consideration the higher 
extinction (A# = 1.6 A^), the range of extinctions explored 
(see below and §[5}, and that the 1.644pm filter width is 25% 
narrower than the 2.14/mi filter. The (2,1) S(l)/(1,0) S(l) ra- 
tios are calculated with the continuum subtracted. The error 
propagation results in very large errors at locations where the 
continuum forms a large fraction of the H2 emission. 

From the KSPEC data, we can extract several pieces of in- 
formation. We calculate the differential extinction between two 
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transition lines (in magnitudes) using 



A^2.512 1ogf F '^ 2A2 ) 



(1) 



where A\ and A 2 are the transition wavelengths, Fi and F2 are 
the relative fluxes, gi and g 2 are the upper level degeneracies, 
and Ai and A2 are the spont aneous electric qu adrupole transi- 
tion probabilities taken from lWolniewicz et alJ {l998). 

When both lines originate from the same upper level and a 
K-band spectral index of 1.7 is adopted, the absolute extinction 
is given by 



A K = 



(2) 



Three pairs of H2 v=(l,0) lines from the KSPEC data were used 
in order to determine A K : Q(2)/S(0), Q(3)/S(l) and Q(4)/S(2). 
The results are presented in Tabled 

Table 3. A% values (in magnitudes) which have been deter- 
mined from the KSPEC relative fluxes in TableQ] The flux er- 
rors have been propagated and yield a realistic Ak error esti- 
mate. 





Q(2)/S(0) 


Q(3)/S(l) 


Q(4)/S(2) 


A K east 


2.7 ± 1.5 


1.5 ± 0.5 


1.8 ±0.5 


A K west 


3.9 ± 2.3 


2.4 ± 0.5 


3.2 ± 2.0 



Immediately evident is that the extinction is higher in the 
western knots although the extinction determination in the west 
is rather tentative due to the observing conditions. A higher ex- 
tinction than in the east is plausible since the western lobe is 
redshift ed and the embedd ed cloud (observed in H 13 CO + by 
iGueth & Guilloteaul dl999t) ~) extends predominantly in this di- 
rection. In order to determine the average value of A K each line 
pair was assigned with the relative line strength as a statistical 
weight. The statistically weighted average Ak values for the 
eastern and western components are 1.8 and 2.9 magnitudes 
respectively. These values fall within the error limits and are 
used to deredden luminosities in §|4] 

In order to interpret the data, we employ the Column 
Density Ratio (CDR) method. We determine the column of gas, 
Nj, in the upper energy level, k Tj, necessary to produce each 
line. We then divide these values by the columns predicted from 
a gas at 2,000 K in local thermodynamic equilibrium with an 
ortho to para ratio of 3. Normalising to the (1,0) S(l) line, we 
then plot this as the CDR, as displayed in Fig.|4] We immedi- 
ately see that the CDR is not constant but a function of excita- 
tion temperature, and thus is not consistent with emission from 
a uniform temperature region or from a single planar shock. 
Furthermore, no significant deviations from a single curve are 
identified (apart from that derived from the (1,0) S(5) which is 
not confirmed by the (1,0) Q(7) from the same Tj). Hence, the 
ortho to para ratio is consistent with the value of three, usually 
associated with H2 shocks. 
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Fig. 4. Column density ratio diagrams for HH211, produced 
from vertical slits running through locations of peak emis- 
sion. The extinction has been adjusted to minimise the differ- 
ences in 1-0 S-branch and Q-branch lines originating from the 
same upper energy level. H2 (1,0) transitions are represented 
by squares, (2,1) transitions by crosses, and (3,2), (3,1) and 
(3,0) transitions by triangles. The faint squares represent (1,0) 
Q branch measurements. 



The extinction determined in the east location also con- 
strains the H-band lines. A significantly lower extinction would 
raise the CDRs for these ( 1 ,0) H 2 above that of the (2,1) K-band 
lines which would have implied non-LTE low density condi- 
tions. As it stands, the fact that a single curve is predicted sug- 
gests a density sufficiently high to ensure LTE. Given a high 
fraction of hydrogen atoms, the lower vibrational levels of hy- 
drogen molecules reach LTE at densities above ~ 10 4 cm" 3 . 



4. Analysis 

4.1. Modeling the bow shocks 

The structure of the western outflow, shown in detail in Fig. [5] 
raises an ideal interpretation scenario. The excited H2 can 
be found in several disti nct knots along an outflow axis. 
iMcCaurfirean et al 1 ill994l) suggested that this well organised 
appearance might prove particularly amenable to modeling. A 
series of bow-shocks is propagating along the outflow. They 
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Fig. 5. The HH211 western outflow is shown here in H2 
(1,0) S(l), H2 (2,1) S(l) and [Fell] emission lines. The images 
are not continuum subtracted as the bow shocks are relatively 
pure. Greyscale bars represent flux levels per arcsec -2 in units 
of 10~ 18 W irT 2 . 



gradually exhaust their momentum and slow down as they 
plough through either ambient gas or the material in the wakes 
of upstream bow shocks. They encounter less dense gas to- 
w ards the edge of the cloud. The submillimeter observations 
of lChandler&Richerl feOOOh show that the HH 211-mm en- 
velope density decreases with distance from the source. The 
azimuthally averaged density structure is well fitted by a sin- 



gle power-law, p < 



-1.5 



, out to 0. 1 pc (the projected distance 



of the outer knots of the outflow from the source). However, 
the density profile is extended in a direction perpendicular to 
the outflow and the outer western knots are located outside this 

dense core region. 

The bow-shock models described by ISmith etafl J2003t) 
have been successfully employed to interpret the se- 
quen ce of bow shape d features along the HH240 outflow 
JO'Connell et alJl2004f) . Here we apply the same model which 
serves as a useful interpretive device. We wish to determine 
if a systematic change of one or more parameters results in a 
close match with the observed bow structures and thus to anal- 



Table 4. Model parameters derived to fit the bow images with 
C-type shocks. 



Parameter 


bow-iie 


bow-fee 


bow-a 


Size, L bm . (cm) 


— 

1.0 x 10 16 


— 

1.0 x 10' 6 


— 

1.0 x 10' 6 


H Density, n (cm -3 ) 


8.0 x 10 3 


4.0 x 10 3 


3.0 x 10 3 


Molecular Fraction 


0.2 


0.2 


0.2 


Alfven Speed, va (kms -1 ) 


4 


4 


4 


Magnetic Field QiG) 


193 


137 


118 


Ion Fraction, x 


1.0 x 10~ 5 


3.0 x 10~ 5 


5.5 x 10- 5 


Bow Velocity, v hm . (kms -1 ) 


55 


40 


29 


Angle to l.o. s. 


100° 


100° 


100° 


s Parameter 


2.10 


1.90 


1.75 


Field angle, /u 


60° 


60° 


60° 



yse the outflow's changing environment. We have applied the 
bow-shock model to the three leading structures: de, be and a. 

Jump-type (J-type) bow shocks can explain many of the 
features associated with high excitation emission regions. 
They cause rapid heating and disso ciate H? mo lecules for 
shock speeds greater than 24 km s~' (Kw anll 19771) . However, 
Continuous-type (C-type) bow shocks have proved extremely 
successful at interpreting most of the observed structures. The 
measured low fraction of ions in molecular clouds is con- 
sistent with their application. The magnetic field cushioning 
means that less energy goes into molecule dissociation; they 
can explain the high infrared fluxes which are observed in bow 
shocks. However, their observed velocities ofte n exceed th e H2 
dissociation speed of ~ 40-50 km s 1 dSmithet alJl99lb giv- 
ing rise to a double-zone bow shock composed of (1) a curved 
J-type dissociative cap (responsible for atomic emission) and 
(2) C-type wings (where H2 emission is radiated). 

Our C-type bow shock model consists of a three dimen- 
sional curved surface described by 



Z/L how = (l/s)(R/L bow Y, 



(3) 



where Z and R are cylindrical coordinates, L\, ow characterises 
the bow size and s defines the shape or sharpness of the bow. 
The curved shock surface is then divided into a very large num- 
ber of steady-state planar shocks for which the detailed physics 
and chemistry are computed (Draine 1980; Smith et al. 1991). 
Each mini-shock propagates at a different velocity depending 
on the angle between the shock surface normal and the direc- 
tion of motion of the bow. The temperature reached, and hence 
the excitation conditions, depends on the shock velocity, den- 
sity, magnetic field strength, etc., as well as the magnetic field 
direction (see Table0J. 

The sy stematic method of ex ploration of parameter space 
is given in lO'Connelletail ^2004). In summary, the bow lu- 
minosities provide constraints on the density and bow speed. 
The location of the emission in the flanks or apex also con- 
strains the bow speed. In addition, the ion fraction constrains 
the transverse bow thickness as well as the bow speed. The 
magnetic field strength strongly influences the extent of the 
wing emission and the atomic fraction affects the line ratios. 
As demonstrated by previous modelling, the bow speed is lim- 
ited to within ~ 15%, the orientation to within 10° and the den- 
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Table 5. Observed and predicted bow shock luminosities and 
(2,1) S(l) / (1,0) S(l) flux ratios. Luminosities are expressed in 
units of Lq and a distance of 315 pc is adopted. Luminosities 
have been dereddened using K-band and H-band extinctions 
of 2.9 and 4.5 m agnitudes (using A H = 1.56 x A K from 
iRieke & Lebofskvl i 1985b ) although the errors in these values 
are relatively large, see table[3] respectively. Note that the [Fell] 
luminosities have been predicted using a J-type shock model. 



Line 


Observed* 


Dereddened 


C-type 








Model 


bow-de 


H 2 (1,0) S(l) 


4.6 x 10~ 4 


6.5 x 10~ 3 


6.8 x 10- 3 


H 2 (2,1) S(l) 


4.3 x 10~ 5 


6.1 x 10~ 4 


7.3 x 10~ 4 


H 2 (3,2) S(3) 


9.8 x 10- 6 


1.4 x 10~ 4 


1.6 x 10- 4 


[Fell] 4 D 7/2 - 4 F 9/2 


3.6 x 10- 5 


2.2 x 10~ 3 


2.3 x 10~ 3 


2/1 ratio 


0.09(0.02) 


0.09(0.02) 


0.11 


bow-bc 


H 2 (1,0)S(1) 


9.7 x 10~ 5 


1.4 x 10~ 3 


1.4 x 10~ 3 


H 2 (2,1) S(l) 


6.85 x 10- 6 


9.7 x 10- 5 


1.5 x 10- 4 


H 2 (3,2) S(3) 


no det.* 




2.7 x 10- 5 


[Fell] 4 D 7/2 - 4 F 9/2 


no det.* 




1.2 x 10~ 4 


2/1 ratio 


0.07(0.02) 


0.07 (0.02) 


0.12 


bow-a 


H 2 (1,0) S(l) 


2.2 x 10~ 5 


3.1 x 10~ 4 


3.2 x 10~ 4 


H 2 (2,1) S(l) 


1.7 x 10~ 6 


2.4 x 10~ 5 


2.5 x 10~ 5 


H 2 (3,2) S(3) 


no det.* 




2.6 x 10~ 6 


[Fell] 4 D 7/2 - 4 F 9/2 


no det.* 




1.2 x 10~ 6 


2/1 ratio 


0.08 (0.04) 


0.08 (0.02) 


0.08 



'The 3cr detection limits for the observed knots over an 8"circular 
aperture are: 2.6 x 10~ 5 for (1,0) S(l); 1.8 x 10~ 6 for (2,1) S(l); 6.4 x 
10~ 6 for (3,2) S(3); and 3.1 x 10~ 5 for the [Fell] image. 



sity, magnetic field, atomic fraction and ion fraction to within 
a factor of two. 

The observed H2 (1,0) S(l) luminosity for each bow pro- 
vides us with the strongest density and velocity constraint. For 
a given bow size, L b0 w, the line emission is directly propor- 
tional to the mass density (= 2.32 x 10~ 24 n c ) and (v bo w) 3 - The 
observed bow luminosities have been dereddened using a K- 
band extinction of = 2.9 mag (see §[3}. Table [5] lists the 
observed dereddened luminosities for each knot along with the 
predicted model luminosities. 

The parameters selected in order to model each bow are 
given in Table@] A constant Alfven speed is maintained so that 
the magnetic field varies with density (B -va (47rp) 5 ). From the 
H2 radial velocity structure described in LSalasetalJll2003l) and 
that the outflow H2 knots have an average velocity projected 
onto the plane of the sky of ~ 50 km s (such a velocity is con- 
sistent with our modelling results) we have estimated an incli- 
nation angle to the plane of the sky of between 5° and 10°. The 
bows have been modelled propagating at this angle, i.e. 100° 
to the line of sight as the western outflow is redshifted. At this 
angle the model images closely resemble the observations sug- 



gesting that the bow speed is similar to the shock speed, i.e the 
bows are propagationg in a medium which is relatively at rest. 

The field angle, ju, is the angle be tween the bow direction 
of motion and the magnetic field (see lO'Connell et all (|2004) 
for a detailed description of the geometry). A misalignment of 
these directions results in asymmetric bow wings such as are 
observed for HH211. In our case the best results were found 
for// -60° ±10°. 

Fig. [6] displays the simulated bow image to compare to 
bow-de in (1,0) S(l) and (2,1) S(l) ro-vibrational transition 
lines of H2, as well as the [Fell] 4 D 7 /2 - 4 F9/2 transition line. 
The shock speed is 55 km s _1 and the pre-shock density is 8 
x 10 3 ctrT 3 . The bow size, L DOW , has been chosen in order to 
match the distance between the upper and lower wings, in this 
case ~5". 

The [Fell] emission is generated by a J-type dissocia- 
tive bow and is restricted to a compact but elongated zone 
towards the bow apex where the highest temperatures are 
reached. However, the observed [Fell] emission is restricted to 
a single compact condensation, unlike the model distribution. 
According to this model the total cooling in the (1,0) S(l) line 
is 6.8xl0~ 3 Lq, 1.4x10~' Lq in all H 2 rotational and vibra- 
tional lines and 2.4xl0 _1 L0 in all atomic and molecular lines. 

Bow-fee is observed in both the (1,0) S(l) and (2,1) S(l) 
lines. Fig.0displays the model generated images which closely 
resemble the observations. The distance between the upper and 
lower wings is ~3". This knot is modeled with a reduced bow 
speed of 40 km s which is propagating into a lower density 
medium (n = 4.0 x 10 3 cm -3 . No [Fell] emission is observed, 
consistent with the fact that its predicted luminosity lies below 
the detection threshold. According to this model the total cool- 
ing in the (1,0) S(l) line is 1.4xl0~ 3 L , 3.6xl0~ 2 L Q in all H 2 
rotational and vibrational lines and 5.4xl0~ 2 L0 in all atomic 
and molecular lines. 

Bow-a appears in the (1,0) S(l) line as a compact knot of 
emission with a (2,1) S(l) / (1,0) S(l) ratio of 0.08 + 0.02. A 
bow speed of 29 km s results in emission restricted to the 
bow apex. The calculated line cooling is 3.2 x 10 -4 Lq in the 
(1,0) S(l) line, 8.8 x 10 3 L in all H 2 lines and 1.3 x 10~ 2 L Q 
total line cooling. 

We conclude that the observed structures are recreated in 
the model by solely altering the density, bow speed and ion 
fraction while all other parameters remain fixed. The extent of 
the bow is influenced principally by the bow speed as it deter- 
mines the location of the H2 emission. A slower bow is char- 
acterised by emission concentrated closer to the bow front. For 
this reason we have maintained a constant L bow for all the mod- 
els. The bow speed systematically decreases as an otherwise 
similar bow ploughs into a material of decreasing density. An 
increase in the ion fraction is expected in less dense regions 
where cosmic rays and UV photons can more easily penetrate 
the gas. 

The driving power of a bow is converted into heat at a the- 
oretical rate P - lpv\ L? 7r/2 or 

~~ bow bow ' 

P = 13/rC - V Vbmt ) 3 ( Lbow ) 2 L,., 
H8.0x 10 3 cm- 3 /\55kms" 1 /\1.0x 10 16 cm/ °' 

(4) 
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Fig. 6. A C-type bow shock model for bow-de shown in H2 
(1,0) S(l) and (2,1) S(l) excitation lines. Adopting a source 
distance of 315 pc gives a pixel scale along the x and y axes of 
1 pixel = 0".18. The bow direction of motion is inclined to the 
plane of the sky by an angle of 10°, away from the observer, 
i.e. 100° to the angle of sight. The bow is moving at 55 km 
s _1 relative to the ambient medium of H density,n, = 8 x 10 3 
cm" 3 . All flux levels are indicated normalised to the maximum 
H2 (1,0) S(l) line flux. Hot dissociative shocks are necessary 
to induce emission from [Fell]. Here we have employed a J- 
shock model (same parameters) to simulate the flux distribution 
at 1.644 yum. 
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Fig. 7. C-type bow shock model for bow-frc which reproduces 
the observed H 2 (1,0) S(l) and (2,1) S(l) flux distribution. 
Adopting a source distance of 315 pc gives a pixel scale of 
1 pixel = 0".18. Here we take a lower bow speed of 40 km 
s _1 and a pre- shock H density of 4 x 10 3 cm -3 to match the 
observed luminosity. Other parameters are given in table |4] 
The grey scale flux indications are normalised to the maximum 
(1,0) S(l) level forbow-de in Fig.© 



where £ is a non-dimensional factor related to the aerodynam- 
ical drag and is of order unity, p is the mass density equal to 
2.32 x 10~ 24 n (n is the hydrogen nucleon density). The derived 
numbers are thus consistent with expectations. 

4.2. The outflow continuum emission and excitation 

We have also divided the fluxes in the narrow -band images con- 
taining the (1,0) S(l) flux and the (2,1) S(l) flux to produce the 
ratio image displayed in Fig. [5] In order to produce this map 
both images were first smoothed with a Gaussian FWHM 
()".6 and values lying below the noise level were excluded be- 
fore dividing. Note that the images were not continuum sub- 
tracted, due to the poor quality of the continuum image at 
2.14yum, so Fig.|5]reveals two distinct outflow regions: (1) areas 
where the continuum emission is relatively strong possess a ra- 
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Fig. 9. A map showing the ratio of our 2. 122yum and 2.248yum images. The two images have not been continuum-subtracted (due 
to the very low S/N of the continuum image), so the map only reveals the H2 (2,1) S(l) /(1,0) S(l) line ratio towards the ends of 
the outflow, where the continuum is weak; the darker regions near the center of the image are where the emission is dominated by 
the continuum. Logarithmically increasing contour levels are at 0.08, 0.1 1, 0.16, 0.22, 0.32, 0.45 (black) and 0.63, 0.89 (white). 



tio above 0.3 and approach 1.0 where only continuum emission 
is present and (2) parts of the outflow not effected by continuum 
emission where the H2 excitation (~ 0.1) can be studied. 

The outflow regions dominated by continuum emission 
are located along the e dges of the SiO jet imaged by 
IChandler & Richenll200ll) . This supports the idea that the con- 
tinuum arises through scattered light from the protostar. The 
light escapes along a jet excavated cavity but not along the high 
density jet itself. It is cut off at 17"(8.0 x 10 16 cm) along the 
western outflow where it possibly encoun ters a high density 
1 Mq filament (Gueth & Guil loteaull 1 9991) . It is also possible 
that the filament lies in front of the flow. This projection effect 
would explain the higher extinction measured in the western 
outfl ow, contrary to the submillimeter dust emission maps of 
llChandler & Richeifc oOO) which do not reveal a strong asym- 
metrical density distribution between the eastern and western 
outflows. Continuum light encounters less hindrance along the 
eastern outflow where it terminates alongside shock excited H2 
emission at knots i and j, 45"(2.1 x 10 17 cm) from its pro- 
tostellar origin. Here continuum emission is seen at 2.14/im 
(Eisl offel etal.l 2003) but not in our ratio image as the H2 emis- 
sion is relatively strong here. The origin and implications of the 
continuum emission will be discussed in §|5] 



The excitation ratio along the western outflow,- 0.1, is 
typical of outflows seen in col lisionally excited emission 
(Black & Dal garnolll976l IShull & Beckwithlll982h . There is 
an increase in the excitation ratio towards the leading edge 
of bow-iie where higher shock velocities and temperatures are 
reached and the H2 line emission reaches its maximum value. 
Similar conditio ns wer e found for the HH 240 bow shocks 
(O'Conn ell et al.l 12004) and they are well explained through 
the bow shock interpretation. However, a strong deviation from 
this picture is found in the [Fe II] image: The localised [Fell] 
emission is coincident with a region of higher ratio, (labelled 
[FeII]-3) and, puzzlingly, not at the expected H2 dissociated 
bow apex as seen in the model generated image, Fig. [6] 

5. Discussion 

Which mechanisms give rise to the series of near-infrared bow 
shocks? Our analysis suggests a combination of two processes: 
(1) The bow shocks are generated by a series of similar out- 
flow accretion/ejection events. Using the model velocities the 
time lapses between the bows (and therefore between outflow 
events) are ~ 425 and ~ 290 years. These numbers are consis- 
tent with the detection of three bows given the dynamical age 
of ~ 1000 years. Fluctuations in jet activity probably manifest 
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Fig. 8. Faint (1,0) S(l) emission is detected at bow-a. A bow 
propagating at 29 km s _1 generates the correct (1,0) S(l) lu- 
minosity with luminosities from the other lines below the level 
of detectability, see table |5] Adopting a source distance of 315 
pc gives a pixel scale of 1 pixel = 0".18. The grey scale flux 
indications are normalised to the maximum (1,0) S(l) level for 
bow-de in Fig. [6] 



themselves as bow shocks at the jet cloud impact region which 
is Knot-/ where the high speed CO jet terminates and H2 shock 
heating is initialised. At Knot-ftw o H? (1,0) S ( l) velo city com- 
ponents have been observed by Sal as et alJ (J2003). This re- 
verse/transmitted shock pair seems to indicate the critical zone 
of jet impact where H2 bow shocks are born. After formation, 
the bows propagate away from the protostar towards the cloud 
edge and through a changing environment; they become less 
luminous as the density decreases and they loose their momen- 
tum. Bow-a represents the final stages in the detectable life of 
one of these bow shocks. In support of this process, we note that 
the molecular jet demonstrates an apparent acceleration along 
its len gth of value ~ 5 x 10~ 3 km s -1 ALT 1 (Chandler & RicheJ 
2001). Given ballistic motions, this implies that the entire jet 
now observed was ejected within a relatively short period of 
time just 100-150 years ago. 

(2) The bow shocks become illuminated within regions 
where the outflow impacts on denser clumps of gas. This idea 
is supported by the implied high K-band extinctions of 2.9 and 
1.8 magnitudes for knots de and i. These values do not repre- 
sent the entire outflow and cannot be used to infer the dered- 
dened luminosity for the whole outflow. Additionally, strong 
continuum emission is detected in the eastern outflow at knots 
i and j where the extinction is high (A K ~ 1.8). If this is indeed 
scattered light from the protostar then it has channeled through 
to where the enhanced density has resulted in scattering. The 
passage of a C-shock will increase the post-shock density with 
a compression ratio of y2 times the magnet ic Mach num - 
ber (shock speed divided by the Alfven speed) (Snit zerll978l) . 
Therefore, the continuum emission is likely to be seen along- 
side the shock excited H2 emission as is the case for HH 211. 



iMcCaughrean et al suggested an average A% of 1.2 

magnitudes which is lower than the values that we derived in 
§ 3 but is consistent to within our error limits. Adopting this as 
the average extinction over the entire outflow we find that the 
total H 2 (1,0) S(l) luminosity of the entire outflow is 0.009 L©. 
According to the predictions of our bow shock models, ~ 4% of 
the total H2 line emission (2.7% of total line cooling) is emit- 
ted in the (1,0) S(l) which gives an intrinsic H2 luminosity for 
HH 211 of ~0.23 L.0 and a total luminosity, L ra£ /, resulting 
from all line emission of ~0.34 L.0. If all the mechanical en- 
ergy is converted into radiation then L mec i, is equivalent to L ra< i 
and L mech / L bol -10%. 

Provided that the velocity of the swept up gas, which coin- 
cides with the CO outflow, roug hly equals the sh ock velocity 
(i.e. a radiative shock). iDavis & Eisloeffell (11996) have shown 
that L rac { should be roughly equal to the kinetic luminosity Lfa„ 
of the outflowing material as measured through the CO lumi- 
nosity. These criteria are indeed met for HH 211 as L/. m ~ 
0.24 L.0 and the CO and H2 are coincident su ggesting that 
the shock is essenti ally radiative JGueth & Guilloteaul ll999: 
iGiannini et alJEoOll) . Note that the average extinction is re- 
stricted to about 1 magnitude in the K-band in order to meet 
these criteria. 

The driving source HH 21 1-mm has a bolometric tempera- 
ture of 33K which implies a youthful outflow system. However, 
the HH 211 outflow itself does not reveal any characteristic 
signs of its assumed youthfulness besides the small spacial ex- 
tent which implies a dynamical timescale of order 1000 years. 
This outflow age is limited by the density structure of the en- 
vironment, as older bow shocks may simply have disappeared 
into sparse material. For this reason the outflow extent itself 
cannot be used to infer the full duration of outflow activity. 

The high density jet (2-5 x 10 6 cirT 3 ) observable in SiO J 
= 5— > 4 emission implie s a pre-shock density of about 2 x 10 
cm " 3 (iGibb et alj r2004) and thus a maximum jet-to-ambient 
density ratio of ~20. The CO J _ 2— >1 maximum radial ve- 
locity is ~ 40 km s~' JGueth & Guilloteaulll999l) implying a 
jet velocity of 230-460 km s _1 given an inclination angle to 
the plane of the sky of between 5° and 10°. Hydrodynamic 
numerical simulations find that a low density envelope gener- 
ally develops around jets of density ~ lO^cm" 3 and a jet- to- 
ambient density ratio of 10 JSuttner et alJI 19971 IVolker et alJ 
ll999HRosen & Smithl2004l) . Could (pulsed) C-type jets which 
are heavier, denser and more ballistic than those simulated ex- 
cavate low density cocoons through which continuum emission 
might escape along the outflow? To date, such high density and 
high velocity jets involving C-type physics have not been simu- 
lated. The new challenge for theorists is to simulate such MHD 
jets. 

HH 21 1 is of particular interest due to its unusually strong 
continuum emission. The original radiation may have escaped 
from the source along a cavity of low optical depth which must 
have been excavated by powerful jet events. This radiation is 
then scattered when it encounters dense walls or clumps where 
the NIR optical depth along the direction facing the source is 
of order unity. Some of the scattered radiation then exits along 
the line of sight through a foreground of moderate NIR extinc- 
tion. This interpretation thus requires the original radiation to 
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penetrate a distance of order 10 17 cm with a mean density of 
less than 10 5 cirT 3 before encountering walls of thickness of 
order 10 16 cm and density 10 6 cirT 3 . A moderate fraction of 
the scattered radiation then escapes without encountering fur- 
ther dense features along the line of sight. The high density 
and dumpiness of the jet as revealed throug h SiO observations 
JChandler & Richerl200lHNisini et alJ2002h are important fac- 
tors to consider in this interpretation. 

The [Fell] emission at 1.644 /mi originates from an upper 
energy level of 1 1,300 K (compare to 6,953 K for H 2 (1,0) S(l)) 
therefore we expect it to highlight the hotter, high excitation re- 
gions of an outflow. As shown in § |4] the [Fell] emission from 
a typical bow shock should be located towards the front of the 
bow, where H2 is dissociated. Clearly, this is not the case for 
HH 21 1 (bow-ofe). The most likely explanation is that the front 
of the bow shock is traversing a low density region. Material 
has been swept out by outflow activity and the bow shocks be- 
come luminous only where they interact with the wall of this 
hollowed out cavity. In this way it is possible see [Fe II] emis- 
sion coincident with bright H2 emission in the 'shoulders' of 
the bow shock. Such an explanation is also supported by the 
distinct possibility that continuum emission from the source is 
escaping along an outflow excavated cavity. 

6. Conclusions 

We have studied the HH211 protostellar outflow in the near- 
infrared regime through high resolution imaging and spec- 
troscopy. Images in the (1,0) S(l) and (2,1) S(l) ro-vibrational 
transitions of H2 have been analysed in order to study the 
excitation throughout the outflow. A narrow-band image at 
1 .644/mi was presented where the outflow is clearly seen and 
several confined condensations of emission are identified. In 
addition, we have presented K-band spectroscopic fluxes for 
two separate prominent locations from which the extinction 
and excitation conditions have been investigated. We have suc- 
cessfully modeled the series of bow-shocks in the western out- 
flow as curved 3-dimensional shock fronts with steady state 
C-type physics. 

Our findings have lead to several conclusions about the na- 
ture of the outflow: 

- C-type bow shocks propagate along the western flow. 
Model fitting has constrained several parameters including 
the density, bow velocity, ion fraction, intrinsic bow shape 
and magnetic field strength and direction. 

- The bow shocks are passing through dense clumps where 
their luminosity is accentuated. High values of extinction 
are measured in these regions. 

- Bows de, be and a are plausibly modeled as a series of ini- 
tially identical bow shocks propagating through a medium 
of decreasing density. The bows slow down and disappear 
as they approach the cloud edge. 

- A misalignment of the magnetic field and outflow direc- 
tions can account for the observed bow shock asymmetries. 

- The most likely source of the protracted continuum emis- 
sion is light from the protostar which evades dense core 
obscuration by escaping through a low density cavity exca- 



vated by the jet, as discussed in § [5] The continuum light 
is scattered when it encounters the denser material aligning 
the jet tunnel and the dense clumps along the outflow. 

- The excitation along the outflow is typical of outflows in 
general. The ortho to para ratio of 3 for molecular hydrogen 
implies shock heating as the source of the near-infrared line 
emission. 

- The [Fell] emission is predicted and detected in isolated 
condensations. These condensations are coincident with 
strong H2 emission. However, the location of the [Fell] 
emission is puzzling; it is not found in the expected bow 
apex region as predicted. This may also be due to the low 
density tunnel through which the bow apex is propagating. 

These findings together with the large volume of previously 
published material is suggesting a global outflow model, as fol- 
lows. Episodic fluctuations in accretion/ejection (of order a few 
hundred years) give rise to a variable jet velocity. The resulting 
shocks manifest themselves as C-type bow shocks at the princi- 
ple jet/ambient medium impact region where they are detected 
outside the dense core where the extinction is lower. The bows 
propagate towards the cloud edge through a changing environ- 
ment. The model suggests that the mean density decreases with 
distance from the core but that the bow shocks brighten where 
they encounter dense clumps. It is feasible that these clumps 
consist of gas swept up by passage of previous bows driven by 
the alternating outflow power. 

It is clear that in-depth studies of a wide range of protostel- 
lar outflows will yield valuable insight into how the cloud en- 
vironment sculpts the outflow and how much the environment 
itself has been influenced by the star forming process. 
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